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Abstract: Some regions of parameter space of the minimal supersymmetric standard 
model (MSSM) with high scale supersymmetry breaking have extreme sensitivity of elec- 
troweak symmetry breaking (EWSB) to the top quark mass through renormalisation group 
evolution effects. This leads to uncertainties in the predictions which need to be taken into 
account in the interpretation of searches for supersymmetric particles in these regions. As 
an example, we provide estimates of the current uncertainties on the position in parameter 
space of the region which does not break electroweak symmetry in the constrained MSSM 
(CMSSM). The position of the boundary of EWSB can vary by up to 2 TeV in m,Q due 
to the uncertainties coming from the current measurement errors on the top quark mass 
and from higher order corrections. In this dangerous region, for fixed CMSSM parameters 
the neutralino lightest supersymmetric particle mass has an associated large uncertainty 
of order 100%. These uncertainties therefore have a profound effect on the interpretation 
of LHC supersymmetric particle searches in terms of the CMSSM. We also show how to 
ameliorate poor convergence of the iterative numerical algorithm that calculates the MSSM 
spectrum near the boundary of EWSB. 
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1 Introduction 

Searches for supersymmetric particles are one of the highest priorities for the Large Hadron 
Collider, and the subject of many recent experimental and theoretical publications. The 
data taken up until 2011 at a centre of mass energy of 7 TeV has now been extensively 
analysed. New results have begun to appear from 8 TeV data, and results from the complete 
8 TeV data set are expected to appear in 2013. The most stringent limits on the masses of 
supersymmetric squarks and gluinos come from generic searches using the missing energy 
signature, with a varying number of high transverse momentum jets [1, 2]. The searches 
in Ref. [1] exclude gluino and squark masses below 860 GeV and 1320 GeV respectively in 
simplified scenarios where only direct decays of these sparticles are allowed and R— parity is 
assumed. In the constrained MSSM (CMSSM) framework [3-7], where the sparticle masses 
are fixed at a high gauge unification scale Mx ~ 2x 10 16 GeV, the searches exclude universal 
gaugino masses Mu2 < 640 GeV at low values of universal scalar mass mo, and M 1 / 2 < 
300 GeV at high mo, approximately independently of the ratio of the two Higgs vacuum 
expectation values tan/3 = V2/V1 and the trilinear supersymmetry (SUSY) breaking scalar 
coupling Aq. In the CMSSM the limits approach the region where electroweak symmetry 
breaking (EWSB) does not occur, and the search limits are generally not plotted beyond 
the expected boundary of this region. The same can be seen in the most recent conference 
interpretations of so far unpublished data from these experiments. 

When the LHC runs at higher energy (13 or 14 TeV, starting in 2014), the searches will 
be sensitive to far higher values of mo and -Mi/ 2 - It is therefore important to understand 
the behaviour of the theoretical predictions as the boundary of EWSB is approached. The 
uncertainty on the location of the boundary needs to be taken into account, as well as the 
uncertainties on the predicted sparticle masses. These effects can significantly alter the 
region of validity of the exclusion limits in interpretations of the data in terms of high- 
scale SUSY breaking models (although they of course do not alter the data themselves). 
In particular, if one is attempting to rule out the CMSSM (i.e. a bona fide top-down 
analysis), the question of the uncertainty of the position of the electroweak symmetry 
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breaking boundary is a physical one: in this case, we must fix the CMSSM parameters 
and ask if the point is ruled out or not. This is also in principle the case for other top 
down high-scale models of supersymmetry breaking where the various SUSY breaking 
terms are set by a small number of parameters. If one is performing a 'simplified model' 
or 'phenomenological MSSM' approach, one has enough free parameters to make sure 
that EWSB is broken correctly, and the problem should not arise. In practice, there is 
much interest from the experiments in ruling out the CMSSM in particular because it is a 
familiar and a priori reasonable model, and because it serves as a well-studied benchmark 
for experimental progress in model space. 

For a given point in MSSM parameter space, the computer program SDFTSUSY [8] 
is currently used by ATLAS and CMS to calculate the sparticle spectrum and MSSM 
couplings, and to determine whether or not the point performs EWSB successfully. We 
shall therefore also utilise SOFTSUSY to quantify uncertainties due to the position of EWSB. 
Near the boundary of successful EWSB, SOFTSUSY reports non-convergence of the numerical 
algorithm, which manifests as a small parameter region on the experimental plots where 
the interpretation of an experimental search in terms of a high-scale SUSY breaking model 
cannot be trusted. Refs. [9] and [10] have also recently noticed that at high values of third 
family scalar masses, sparticle spectrum calculators (including SOFTSUSY) give inaccurate 
results: the sparticle mass calculation is particularly difficult to solve in this tricky region 
of parameter space. The region of parameter MSSM space with large scalar masses has 
received considerable attention in the literature, and is called the focus point, or hyperbolic 
branch [11—13]. The focus point has received renewed attention recently because if the 
SUSY breaking scalar trilinear coupling involving Higgs and stops At is large, then a 125 
GeV lightest CP even Higgs boson (as implied by 2012 LHC Higgs search data [14, 15]) 
may result [16-18]. 

1.1 EWSB in the MSSM 

We write the MSSM superpotential in terms of chiral superfields, suppressing gauge indices 
and neglecting all Yukawa couplings save those of the third family, 

W = nHzHx + h t Q 3 H 2 u 3 + h b Q 3 H 2 d 3 + h T L 3 H\e 3 , (1.1) 

where the chiral superfields of the MSSM have the following Gsm = SU(3) C x SU(2)l x 
U(1)y quantum numbers 

Li: (1,2,-1), e i: (1,1,1), Qi : (3,2, i), Ui :(3,l,~), 

di-. (3,1,1), H x :(l,2,~), H 2 :(l,2,±), (1.2) 

and i £ {1, 2, 3} is a family index. We shall write superpartners of SM fields with a tilde. 

Minimising the MSSM Higgs potential with respect to the electrically neutral compo- 
nents of the Higgs vacuum expectation values, one obtains the well-known tree-level result 
for the Higgs mass parameter fj, in the modified dimensional reduction scheme (DR) 

9 tan 2/3 r 2 2 •, Ml 
H = — - — [m H2 tan/3 - m Hl cot/3] — . (1.3) 
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In order to reduce 1 missing higher order corrections, all quantities in Eq. 1.3 are understood 
to be evaluated at a DR renormalisation scale Q = Msusy , where Msusy is the geometric 
mean of the two stop masses, tan/3 = (H®) / (H®) is the ratio of the two MSSM Higgs 
vacuum expectation values and mHi 2 are the soft SUSY breaking DR mass terms of the 
Higgs doublets. If m 2 Hi and rn 2 H2 and tan 2/3 are such that fi 2 > results from Eq. 1.3, the 
model point may break electroweak symmetry successfully. On the other hand, if fi 2 < 0, 
electroweak symmetry is not broken successfully and the model point is ruled out. We shall 
here study current uncertainties in the position in MSSM parameter space of the fi = 
contour that separates the region of successful EWSB from the region of no EWSB. We 
shall also study the extent to which these matter in terms of searches for super symmetric 
particles. 

The m/jj. terms on the right hand side of Eq. 1.3 are connected to the squark and 
gluino masses in models of SUSY breaking through the boundary conditions on SUSY 
breaking terms, imposed at a (typically high) scale Mx and through the renormalisation 
group equations (RGEs), which at one-loop order are [19] 



„dmi, 
16tt 2 — ^ = 6 



dt 



16** 



(m 2 H2 + m| 3 + m? 3 + A 2 )h 2 ] - 6g 2 2 M 2 - \g 2 M 2 + -gj (m% 2 - m 2 Hl 

,2 2 r.2 , 2 , 2" 



dt 



Tr[mJ - m\ - 2m£ + m| + mj]) , (1.4) 

(m 2 Hl + m| 3 + m| + A 2 b )h 2 ] - Qg 2 M 2 - ^g 2 M 2 - p\ (mjj 2 - m 2 Hl + 
Tr[m| - m\ - 2m? + m| + m 2 ]) , (1.5) 



where t = InQ. When (mjj 2 + fn 2 Q + m| a + A 2 ) is large, the renormalisation of mji 2 (and 
therefore its value at Msusy), is extremely sensitive to ht, the top Yukawa coupling. h t 
is usually the largest dimensionless coupling of the MSSM and is determined at the weak 
scale by the experimental input for the top quark mass 

h t {M Z ) = . r . , r , 1.6 

v(M z ) sm /3(Mz) 

where v(Mz) ~ 246 GeV is the Higgs vacuum expectation value (VEV) that would be 
extracted for a single Higgs boson (in terms of the two MSSM Higgs VEVs vip, v = 
^vl + v 2 ). The running top mass mt(Mz) is determined from the pole top mass m v ° e by 

m t (M z ) = mf e - Am t (M z ), (1.7) 



where Amt(Mz) are loop corrections to the top quark mass evaluated at the DR renor- 
malisation scale Mz- ht is around 1 for typical MSSM spectra (for instance, at parameter 
point CMSSM10.1.1 [20], h t (M z ) = 0.90). Eqs. 1.4,1.5 are used, along with the other 
RGEs of the MSSM, to evolve 2 along with the other soft SUSY breaking parameters 
to Msusy, where the Higgs potential is minimised and Eq. 1.3 is applied. The region of 
successful EWSB is given by ^{Msusy) 2 > in Eq. 1.3, i.e. (since tan/3 > 1 for models 



1 This prescription at least ensures that the dominant terms do not involve large logarithms. 
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in which the lightest CP even Higgs is heavy enough to satisfy experimental constraints, 
tan 2/3 < 0), 

M 2 

[m 2 H2 {M S usY) tan P - m 2 Hi {M SUSY ) cot 0\ < ^J-. (1.8) 

Eq. 1.8 splits the MSSM parameter space into two regions: one region has successful EWSB 
whereas the other does not. 

In models where we consider the primordial SUSY breaking terms to be set at some 
higher scale M x > Msusy, we shall now show that whether or not this condition is sat- 
isfied depends upon the other soft SUSY breaking parameters as well as m 2 H2 (M x ), since 
m 2 H .{MsusY) is determined by solving Eqs. 1.4,1.5. To leading order in A = log(Mx /Msusy) / (167r 2 ) , 
we solve them to obtain 

m 2 H2 {MsusY) = m 2 H2 (M x ) - A (g \m 2 H2 {M x ) + m| 3 (M x ) + m| 3 (M x )+ 

A 2 (M x )]h 2 (M x )). (1.9) 



mjj^MsusY) = m 2 Hi {M x ) - A (6 [(m^Afx) + m^(M x ) + mf 3 (M x )+ 

A 2 (M x )}h 2 (M x )]), (1.10) 

where we have neglected terms proportional to electroweak gauge couplings g\ 2 in favour 
of those proportional to h 2 or h 2 . This is a good approximation when 6Ah 2 b (M x ) <C 1, 
which is not the case in the CMSSM. However, it is qualitatively instructive: substituting 
Eqs. 1.9,1.10 into Eq. 1.8, we obtain an inequality on the soft SUSY breaking parameters 
evaluated at M x : 



m Ho < TTTT \ 6Ah t ( m % + m i + A t) + 7 7TT o Q + 

2 1 - 6 Ah 2 { 1 \ Qa Ui V tan (5 tan 2/3 

4 - 6Ah 2 (m 2 Q3 + m| + A 2 )] X , (1.11) 



H 

1 



tan 2 /3 

where h t h, m 2 u , rn% , m? , m 2 - , A t and Ah are all evaluated at M x . For large enough 
values of m 2 H2 (M x ), we see that the inequality is not satisfied and EWSB is not broken 
successfully. Qualitatively, this behaviour could also be analytically derived in the CMSSM 
by examining the complicated re-summed logarithm solutions for the soft masses in Ref. [21] 
in the limit that hh is also neglected in the RGEs (although we do not include them here 
because they are not especially illuminating). 

Much emphasis has been placed in the past upon fine-tuning associated with the region 
close to the EWSB bound [12, 13] where m 2 H2 {MsusY) is somewhat insensitive to the value 
of m 2 H2 {M x ) assumed. In Ref. [22], a study of the uncertainties of the position of the 
EWSB boundary in CMSSM parameter space was included, explicitly showing its extreme 
sensitivity to the top quark mass. Since then, however, the experimental precision with 
which the top mass has been measured has been much improved, reducing the concomitant 
uncertainty. The time is ripe for a reanalysis of the uncertainties on EWSB, so that 
the results may be used for future LHC searches. In section 2, we shall find that the 
uncertainties are still large despite the improved precision upon the top mass measurement, 
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and furthermore we shall demonstrate that they impact on the prospects of discovery. In 
particular, such uncertainties should be included in the theoretical error when interpreting 
experimental results in terms of high-scale SUSY breaking models. As a by-product of 
our study, in section 3, we shall improve the numerical evaluation of the supersymmetric 
spectrum in this tricky region of parameter space, as determined by SDFTSUSY3.3.4. In 
section 4, we summarise our results, as well as provide some simple recommendations for 
interpretations of SUSY searches in terms of high-scale SUSY breaking models that, if 
followed, will take the uncertainties inherent in EWSB into account. 

2 Uncertainties in Electroweak Symmetry Breaking 

We now quantify the uncertainty inherent in the position of the boundary of successful 
EWSB, and its effect on interpretations of searches for supersymmetric particles in terms 
of high-scale SUSY breaking models. This will allow us to suggest some practical steps 
for interpreters of data that allow such theoretical uncertainties to be taken into account. 
We expect many models of SUSY breaking in the MSSM to have regions of parameter 
space where EWSB is not successfully broken. There are however some very constrained 
MSSM models of SUSY breaking where EWSB is always successful (for example, the large 
volume string compactified models [23] , where the ratio of the universal scalar mass to the 
universal gaugino mass is constrained to be 1/V3), or phenomenological models where one 
fixes rrijj independently of the other soft SUSY breaking parameters, discarding the no 
EWSB region. However, for most models, we expect problematic regions of unsuccessful 
EWSB to exist, as in the CMSSM. Bad EWSB regions have been found in minimal anomaly 
mediation, gaugino mediation and general gauge mediation models to name just a few (see 
Refs. [24-26], respectively). Here, we shall focus for definiteness on the familiar example 
of the CMSSM, bearing in mind that very similar effects will be present in most models 
of SUSY breaking. We expect that our qualitative conclusions will hold in other models 
which possess such regions in SUSY breaking parameters which are set at a high scale. 
At the very least, our recommendations of how to quantify such uncertainties should be 
followed in any model which has a region of parameter space which has no EWSB. 

In the CMSSM, the focus point region is close to the region of unsuccessful EWSB and 
shows extreme sensitivity to the top quark mass through RGE effects. In this region of 
CMSSM parameter space, {m 2 H + m 2 A + m| + A 2 ) is large in Eq. 1.4, and rriH 2 (MsusY) 
in Eq. 1.3 is very sensitive to the value of hf. As explained above, ht(Mz) is in turn 
fixed by Eq. 1.6, via the running top mass mt(Mz)- This is affected by the pole mass 
input, which has significant experimental errors (we take rat = 173.5 ± 1 GeV [27]) as 
well as higher order corrections which are not included in the approximation. However, if 
one fixes SUSY breaking parameters at the electroweak scale, even if there are constraints 
among the various parameters and a no-EWSB region exists, the dominant uncertainties 
coming from the top mass (which enter via the RGE) are not present. We note that the 
uncertainties originating from mt are distinct from those that originate from two-loop (or 
higher) corrections to the pole masses of supersymmetric particles [28] or higgs bosons [29], 
which are small (typically a few percent). 
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We include all one-loop MSSM contributions to the top mass, as well as the two- 
loop QCD corrections and two-loop RGE corrections. However, two loop SUSY QCD 
corrections to the top mass from squarks and gluinos are not taken into account via Eq. 1.7 
in SDFTSUSY3.3.4, nor are 3-loop RGE corrections to h t . The full MSSM two-loop 0(a 2 ) 
threshold corrections to mt and three-loop RGE effects were calculated numerically in the 
CMSSM [30, 31] to be at the 2 - 3% level for tan (3 = 50 and A = 0. There, it was noticed 
that in the focus point region, the position of the no EWSB forbidden region moves by 
around 300 GeV in rriQ. In order to estimate the size of the 0(a 2 ) threshold corrections for 
our sets of parameters, we use Eq. (62) of Ref. [32], which gives the two-loop SUSY QCD 
contribution in the limit that the squark and gluino masses are all of equal mass M: 
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where Xt = At — \i cot /3. The squark masses are typically much larger than the gluino 
masses in the focus point region, and so this may be a bad approximation to the actual 
correction. However, we may estimate the size of the theoretical error roughly by estimating 
the size of Eq. 2.1 across parameter space (i.e. assuming that the squark and gluino masses 
are all equal to the average of the dimensionally reduced values of both stop masses and 
the gluino mass). 

In our numerical analysis, we interpolate a grid of mo and M]y 2 values for Aq = and 
m b (m b ) = 4.18 GeV by using S0FTSUSY3.3.4 [8]. Fig. 1 shows the estimated size of the 
higher order correction, as determined in the way described above. Across the parameter 
space, we see that \(Amt)2LSQCD(MsuSY)\ < 0.8 GeV. We may approximately model the 
effect of this theoretical uncertainty by absorbing it into the pole mass input uncertainty. 
If one wanted to do this more precisely, one would renormalise (Amt)2LSQCD from Msusy 
back to Mzi where the running top mass is matched to the pole mass. However, the 
difference between including (Am t )2LSQCD(MsusY) and (Am t )2LSQCD(Mz) is a three- 
loop effect, which we neglect. It is possible that we have underestimated the overall size 
of the theoretical uncertainty: after all, we have not included possible effects from 3-loop 
RGEs [33] or from 2-loop threshold corrections to mt [32]. We shall conservatively add 
the 0.8 GeV theoretical uncertainty linearly to the 2a measurement errors upon the input 
value of mt in order to estimate the total uncertainty in EWSB. 

Having found and quantified the error due to the dominant theoretical uncertainty, we 
may combine it with the mt measurement errors in order to see the effect on the location 
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Figure 1. Estimated size of higher order correction (2 loop SUSY QCD correction) 
{Am t )2LSQCD(MsusY) in the CMSSM for tan/3 = 10 and A n = 0. The correction has a size 
as measured by reference to the colour bar on the right hand side. 

of the EWSB boundary. The result is shown in Fig. 2 for Aq = 0, (a) tan f3 = 10 and (b) 
tan f3 = 50. By varying the experimental uncertainty on mt within 2a of its central value, 
we obtain the yellow (lighter) shaded region showing the uncertainty in the boundary's 
location. The curves marked 'EWSB boundary' are, from left to right, for mt = 171.5 
GeV, mt = 173.5 GeV and m t = 175.5 GeV, respectively. Varying m t by a further 0.8 
GeV in either direction, we obtain the grey bands, providing an estimate of the theoretical 
uncertainty. We also show the expected ha discovery reaches for 4 jets plus (0-3) leptons 
for a 14 TeV LHC and 1 fb _1 of luminosity from Ref. [34]. We see that the location of 
the boundary is uncertain to a huge ~ 2 TeV in mo, and that the discovery regions extend 
well into the region of EWSB uncertainty, irrespective of the size of tan/3, even for only 1 
fb -1 of collected data. We also note that the 2a exclusion contours will be at much larger 
values of mo and Mi/ 2 , and so the problem of the uncertain EWSB boundary will be more 
acute for exclusion than for discovery. This combined with larger 14 TeV data sets implies 
that the EWSB uncertainty is going to be very relevant for SUSY searches in the future. 

In Fig, 3, we show that the size of the uncertainty can make a large difference to 
the production cross-sections of squark (anti-)scjuark a ss and (anti-)squark gluino a S g as 
calculated by PR0SPIN02 [35], using the SUSY Les Houches Accord [36] to transfer infor- 
mation about the MSSM spectrum. The solid lines show the limits of the uncertainty, as 
defined by both the experimental 2a mt uncertainties and the 0.8 GeV theoretical error 
acting in tandem. If, for instance, we assume the CMSSM and constrain Mi/ 2 = 400 
GeV from other measurements, squarks may be highly visible at the LHC (when jointly 
produced with gluinos as shown in Fig. 2b) with large production cross-sections (100 fb), 
or they may require a luminosity upgrade (1 fb or less) for discovery, depending upon the 
locus of the boundary. Of course, this would not affect SUSY discovery because discovery 
would be dominated by production of the lighter sparticles. However, the interpretation of 
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Figure 2. Uncertainties in the position of the EWSB boundary compared to the expected 5a 
ATLAS 14 TeV CMSSM search reach for 1 ftr 1 from Rcf. [34]. The region below the 'EWSB 
boundary' is theoretically inaccessible: the uncertainty in the position of this region due to 2a 
variations in m t is shown by the yellow band. The theoretical uncertainty is marked by the grey 
region and estimated as detailed in the main body of the paper. 



supersymmetric signals in terms of a high-scale SUSY breaking model would be obscure 
unless the uncertainties are taken into account. For example, it may be that the values 
of too and Min inferred from measurements of a sg and other processes appear to be in 
the EWSB region. One would then erroneously rule the CMSSM out. When searches are 
interpreted near the boundary of EWSB therefore, in a high-scale SUSY breaking model 
such as the CMSSM, it is important to generate default SUSY spectra using a m t input 
value that is 2a higher than the central value, so that the high end of the uncertain region 
is included in the search. Increasing mt has a negligible effect on the gluino and first or 
second generation squark masses, so the production cross-sections should not change much. 
However, we show below that it has an effect on the predicted lightest neutralino mass in 
the high-scale model, mt does have a larger effect on stop searches and so here, one would 
need to check the effect's magnitude by taking the variation of the stop masses from mt 
into account. 

On the other hand, we know that /x — > at the EWSB boundary, affecting the mass of 
the lightest supersymmetric particle (LSP), here assumed to be the lightest neutralino. The 
MSSM Lagrangian contains the neutralino mass matrix as — ^"ip 0T Ai J^ow + h.c, where 
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Figure 3. Uncertainties in EWSB correlated with (a) squark-(anti-)squark and (b) squark gluino 
production cross-section at the 14 TeV LHC for tan/3 = 10 and Aq = 0. The two solid lines show 
the total uncertainty on the position of the boundary of EWSB. The total production cross-sections 
(with no cuts applied) are shown by reference to the colour bar on the right hand side. The dotted 
lines show (a) iso-contours of niq and (b) iso-contours of mg, as labelled. 



t/)° = (—ib, —iw 3 , h\, \i2) T and, at tree level, 

(Mi -M z cps w M z sps w \ 

M 2 Mzcpcw -Mzspcw ^ ^ 

-M z cps w M z cpc w -pL 
\ Mzspsw -M z spc w J 

We use s and c for sine and cosine, so that sp = sin/3, cp = cos j3 and swipw) is t ne sine 
(cosine) of the weak mixing angle. The condition [i — > (at the EWSB boundary) results 
in a lightest neutralino mass (equal to the smallest eigenvalue of Eq. 2.2) that is zero at 
tree-level 2 . Including loop corrections to the neutralino mass matrix, the neutralino LSP 
is still very light, typically less than 1 GeV. As the position of the EWSB boundary is so 
uncertain, the position of this negligible neutralino mass occurs for very different values of 
m,Q. Since mo controls the squark mass, the correlation between squark and LSP mass picks 
up large uncertainties. We illustrate the problem by scanning along the Mi/ 2 = 300 GeV 
line for various different values of mt, including an estimate of the theoretical uncertainty 
originating from higher order corrections to mt as described above. From Fig. 4a, we 
see the large uncertainty induced in the lightest neutralino mass along the line. We see 

2 One can easily see this by the observation that adding tan /3 multiplied by the third column to the fourth 
column leads to a column of zeroes in the fj, — ^ limit. The resulting matrix, which has a determinant 
of tan p multiplied by the determinant of the original matrix, must therefore have a zero determinant and 
therefore possesses a zero eigenvalue. 
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Figure 4. Evolution in certain SUSY particle masses approaching the EWSB boundary for tan/3 = 
10, Mi/2 = 300 GeV and A = in the CMSSM. (a) shows uncertainties in the CMSSM predicted 
lightest neutralino mass. The solid line has nit = 173.5 GeV, whereas the upper and lower dashed 
lines have m t = 175.5 GeV and m t = 171.5 GeV, respectively. Additional uncertainty coming from 
the higher order theoretical uncertainty in m t via RGE effects is marked by the grey region, (b) 
shows values of the squark and gluino masses as a function of mo along the M1/2 = 300 GeV line. 



from Fig. 4b that the high m,Q region where the neutralino mass is particularly uncertain 
corresponds to large values of the squark masses > 1.5 TeV. Here, the classic jets plus 
missing transverse momentum signature will come primarily from gg production, where g 
decays through an off-shell squark into two jets and the neutralino LSP. In the centre of 
mass frame of the 700-800 GeV gluino, the gluino must share its rest mass energy between 
the two jets and the lightest neutralino. Having higher neutralino masses therefore has 
the effect of reducing the px of the jets, on average. Since experimental analyses impose 
lower limits on the jet pts in the cuts, the efficiency of the cuts will decrease for higher 
neutralino masses. We therefore expect the theoretical uncertainty in the neutralino mass 
to induce a non-negligible uncertainty in the cut efficiency. 

In Figure 5 we plot data on the acceptance times efficiency (A x e) for the ATLAS 
search based on all hadronic final states with missing energy [1]. These data, for the 7 
TeV sample, are available in the auxiliary information of Ref. [1]. The data correspond 
to a model in which the gluino decays directly to jets and LSP, with the gluino mass 
fixed at 1050 GeV, close to the EWSB boundary. The other SUSY particle masses are set 
high enough for them to decouple. The figures show the fraction of signal events which 
would pass the experimental selections (based largely on the pt of the jets, and the missing 
transverse momentum, and also be successfully reconstructed after detection, as a function 
of the LSP mass. Figures 5a and b show the lowest and highest effective mass selections 
for each jet multiplicity, respectively. The value of A x e drops rapidly at high LSP mass 
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Figure 5. Acceptance times efficiency Ax e for signal events selected by the ATLAS supcrsymmetry 
search [1], using the signature of missing energy with hadronic jets. The different signal regions 
correspond to final state jet multiplicities (n) between 2 and 6, as shown in the key. The signal 
regions are selected with minimum values of m e ff, the scalar sum of the ptS of the n leading jets 
(appropriate to each signal region) with the missing energy. The lowest available m e ff selection is 
shown in (a), and the highest in (b) for various jet multiplicities. 

for all of the signal regions. This is expected, since the fraction of the initial squark or 
gluino energy passing into the observable jets is reduced as the LSP mass rises, and hence 
the fraction passing the pt selections drops. For the lowest multiplicity final states (SRA), 
with the higher effective mass selection, the effect can be seen even at the lowest LSP 
masses, with Axe reducing from 11.7% to 9.8% as the LSP mass rises from to 150 GeV. 
This represents a change in signal event rate of 17%, a non- negligible uncertainty. In the 
lower effective mass samples, Axe rises in the same LSP mass range in some cases (an 
effect already noted in Ref. [37]). For example, for SRA', it increases from 23.5% to 26.1%, 
increasing the signal rate by 11%, and continues to rise for LSP masses up to 300 GeV. 
This may be due to the behaviour of the cut on the ratio of the missing energy to m e s 
applied in this analysis, since the two distributions will evolve in slightly different ways as 
the LSP mass changes. These effects can be expected to continue to occur at higher centre 
of mass energies, with the final drop in A x e moving to higher LSP masses. 

Even if the high-scale SUSY breaking interpretation of an experiment's sparticle search 
is largely independent of the uncertainty induced in the LSP mass, it is possible that other 
sparticle masses in the signal sparticle cascade decay chains are highly dependent and cause 
large theoretical uncertainties in signal cut efficiencies. In particular, chargino masses are 
highly dependent on ^(Msusy)- Also, third generation sparticles may be highly dependent 
on /j,(Msusy) in some (but not all) regions of parameter space. Extra care to quantify the 
theoretical error should be taken in searches for cascade chains involving these sparticles, 
particularly when they are on-shell (since the mass splittings of on-shell decays may have 
a particularly large effect on the kinematics of the decay products). In order to take 
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SUSY radiative corrections to gi(Mz), ht^ T {Mz)- 



Run to Mx- Apply soft SUSY breaking boundary condition. 



convergence 



■Run to Msusy- 



EWSB, solve for \x and m|. 



Run to M z . 



Calculate Higgs and sparticle pole masses. Run to Mz- 



Figure 6. Iterative algorithm used to calculate the SUSY spectrum. The initial step is the 
uppermost one. Msusy is the scale at which the EWSB conditions are imposed, as discussed 
in the text. Mx is the scale at which the high energy SUSY breaking boundary conditions are 
imposed. 



such uncertainties into account, varying mt within its 2a limits, including the size of the 
shift from higher order corrections (estimated here to be 0.8 GeV) in the estimate of the 
uncertainty will suffice, since the various relevant sparticle masses will vary according to 
the concomitant theoretical error. 

3 Convergence 

As one approaches the boundary of EWSB in the parameter space of high-scale supersym- 
metry breaking models, one finds that the various spectrum generators start to struggle to 
produce numerically stable results. This is because of the extreme sensitivity of the RGEs 
and boundary conditions, but is not a problem if one is only fixing the SUSY breaking 
parameters at the weak scale. S0FTSUSY3 . 3 . 4 is essentially solving coupled, non-linear 
ordinary differential equations with two sets of boundary conditions: one set on the gauge 
couplings gi, the Yukawa couplings h t: b, T and the Higgs potential parameters \x and m| at 
the weak scale, and the other set on the soft SUSY breaking parameters at the high scale 
M X - S0FTSUSY3 .3.4 uses an iterative algorithm, as shown in Fig. 3. For S0FTSUSY3 .3.4 
to report that it has solved the problem, a convergence criterion must be satisfied: the 
values of the DR sparticle masses of the current iteration must all be identical to the pre- 
vious iteration's, to some fractional accuracy e. Here, e = 10 -3 . Near the boundary of 
EWSB, h(Msusy) and ms(MsusY) become extremely dependent upon the precise value 
of ht(MsuSY)- However, h t (MsuSY)-, which is dependent on m t (MsusY), varies from it- 
eration to iteration because the loop corrections depend upon the sparticle masses, which 
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Figure 7. Convergence of the numerical solution to EWSB Higgs potential minimisation conditions 
near the boundary of clectroweak symmetry breaking in the CMSSM using S0FTSUSY3 .3.4. Wc 
have set e = 10~ 3 , m t = 173.5 GeV, tan^ = 10, A = 0, M 1/2 = 337.5 GeV and m = 3000 GeV 
for the 'convergent point', or m =3400 GeV for the 'bad convergence' point. 

also vary from iteration to iteration. In parameter space that is sufficiently near to the 
boundary of EWSB, this numerical iteration noise becomes too great for the program to 
find a solution, as illustrated in Fig. 7. We show two points in CMSSM parameter space, 
one with bad convergence, and one with good convergence. We display the predicted value 
of ^(Msusy), which is one of EWSB variables that is most sensitive. We see that after 10 
iterations, the convergent point has found a stable solution for h(Msusy), whereas there is 
no indication that the point with bad convergence would ever converge. In practice, if no 
convergence is achieved after 30 iterations, S0FTSUSY3.3.4 halts the calculation reporting 
a 'no convergence' error. 

If one's scans in parameter space are sufficiently coarse, one may not be aware of this 
problem, as the region of no convergence is fairly thin in the mo-Mi/2 plane. We zoom 
in on the plane in Fig. 8 in order to investigate it. Fig. 8a shows that the yellow (light) 
regions of bad convergence are prevalent close to the boundary of no EWSB. Using Eq. 1.3 
to predict Mz for the numerical solution that issues from SDFTSUSY3.3.4 (whether or not 
convergence is reached), we plot its value in the zoomed region in Fig. 8b. We see that 
throughout this zoomed region, despite the fact that central experimental value of Mz has 
been used at each iteration to constrain the model, the large sensitivity spoils the prediction 
after each iteration and the predicted value is far from the central one of M| eni = 91.1876 
GeV [27]. The diagonal line in Fig. 8a across parameter space shows the line taken in 
Fig. 9 to illustrate how how the numerical spectral predictions of SDFTSUSY3 . 3 . 4 vary 
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Figure 8. Zoom of CMSSM focus point parameter space as predicted by S0FTSUSY3 .3.4 for 
Aq = 0, m t = 173.5 GeV and tan/3 = 10: (a) shows the region of successful EWSB (black), bad 
numerical convergence (yellow) and no EWSB (red). The diagonal line shows the parameter line 
taken to illustrate non-convergence in Fig. 9. (b) displays the predicted value of Mz/GeV, as 
measured by the colour bar to the right hand side. 

while travelling through this problematic region of parameter space. Each line represents 
a different sparticle mass divided by some constant number of GeV, chosen so that the 
ratios are all similar and appear near to each other on the same plot. We see that the 
gluino, stop and pseudo-scalar Higgs masses have a smooth behaviour as the region of no 
EWSB is approached. However, we see that the neutralino LSP mass, which depends upon 
^{MsuSy)-, becomes erratic close to the boundary, corresponding to numerical problems 
and non-convergence of the iteration algorithm. Using the numerical solution at the inferred 
boundary but substituting \x = into the one loop level neutralino mass matrix, we infer 
a mass of 0.1 GeV for the neutralino. The neutralino line would be at m/M ~ 0.01 at 
this point, rather than around 0.6. Some of the erratic non-convergent behaviour can be 
hidden behind the LEP chargino bound of M x ± > 104 GeV [38], but still some noise 
remains unhidden. 

We improve the algorithm in several ways: firstly instead of just 30 iterations, we 
allow a maximum of 100. This helps the program achieve convergence at some parameter 
points. Secondly, we specify that the predicted value of M^ ole must be close enough to the 
empirical central value for convergence to be possible. Thus, the convergence criteria are 
enlarged to include the constraint that \M^ ole /M| eni — 1| < e within the first 10 iterations, 
or \Mz° le /M^ nt — 1| < lOe after the first 10 iterations. Thirdly, there are some sub 
iterations (for example for the solution of Eq. 1.3 for h(Msusy) once loop corrections are 
included) which are required to converge to much higher accuracy, given the sensitivity 
of each iteration to the underlying parameters. Finally, we further smooth the iteration 
of ^{Msusy) if the parameter point is proving difficult. Between iteration 11 and 20, the 
average of the value of [i u from the previous iteration and the one predicted by the EWSB 
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Figure 9. S0FTSUSY3 .3.4 convergence approaching the boundary of EWSB in the CMSSM for 
m t = 173.5 GeV. The plot shows the blue line in Fig. 8 (i.e. Mi/ 2 is varying across the plot as 
well as mo). We have scaled each i th particle mass by a constant mass Mi for clarity. The region 
where electrowcak symmetry is not satisfactorily broken is denoted 'no EWSB'. m ± > 104 GeV, 
in contravention of the LEP bound, in the area marked 'chargino bound'. 

conditions acting on the current iteration /J, ews b is taken. Between iteration 21 and 30, 
S0FTSUSY3 . 3 . 6 takes h(Msusy) = 0.2/i eu;s b + 0.8^ o /d, whereas after iteration 31, it takes 
I^(Msusy) = 0.1/i eu;s b + 0.9fi o id- We find that this helps dampen the oscillations apparent 
in the points that exhibit poor convergence properties such as the 'bad convergence' line 
in Fig. 7. The resulting improvement in convergence is shown in Fig. 10. The region of 
no convergence has significantly diminished, and become much more regular as compared 
to Fig. 8a. Also, in the region of successful EWSB, Mz is predicted to be much closer 
to M| en *, signalling a good solution to the EWSB conditions. These improvements are 
implemented in an improved version of S0FTSUSY (version 3.3.6) [8]. 

4 Summary and Conclusions 

We have studied the uncertainties associated with the boundary of viable parameter space 
in constrained high-scale SUSY mediation models of the MSSM coming from successful 
EWSB. Even with the current accurate measurements of mt, there is an enormous 2 TeV 
variation in the CMSSM boundary's position as measured by the soft SUSY breaking scalar 
mass input tuq. There is also a theoretical uncertainty associated with higher order cor- 
rections to the top mass, equivalent to around 0.8 GeV. Uncertainties on mt translate into 
an uncertainty on ht, which the RGEs are very sensitive to, and to which, ultimately, the 
EWSB region becomes uncertain. There is therefore no concomitant EWSB uncertainty 
in simplified models or other low-scale models such as the phenomenological MSSM. As 
well as providing an uncertainty on the boundary of viable parameter space, there is also, 
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Figure 10. Zoom of CMSSM focus point parameter space with the improved algorithm in SOFTSUSY 
3.3.6 for nit — 173.5 GcV: (a) shows the region of successful EWSB (black), bad numerical 
convergence (yellow) and no EWSB (red), (b) displays the predicted value of Mz/GeV, as measured 
by the colour bar to the right hand side. 



for a given squark and gluino mass (assuming a particular high-scale model), an uncer- 
tainty on the lightest neutralino mass which results in a theoretical uncertainty upon cut 
efficiencies. Currently, these sources of uncertainty are not taken into account when inter- 
preting experimental SUSY searches in terms of particular high-scale models, and should 
be. However, we have shown that they are significant and so must be taken into account 
for any robust interpretation of data in terms of some constrained model of SUSY breaking 
at a high scale. The significance of the uncertainties is likely to increase as the centre of 
mass energy of the collider increases, since the probed region covers more of the uncertain 
region of parameter space. In the CMSSM in particular, such uncertainties should not 
have a significant impact on SUSY discovery, since the uncertainties occur in a regime of 
large scalar masses, so discovery is based upon the production of the other lighter sparticle 
states. But they do have a large impact on whether the CMSSM is a viable model or not 
(as opposed to some competing SUSY breaking mediation model, for instance). 

We now provide recommendations for those wishing to present robust interpretations of 
searches for supersymmetric particles, which one should apply when considering high-scale 
and constrained models of SUSY breaking such as the CMSSM: 

1. Generate default spectra using using a input value that is 2a plus a theoretical 
uncertainty (currently 0.8 GeV) higher than the central value. This ensures that one 
is not erroneously ruling out a region of parameter space for not breaking electroweak 
symmetry correctly. 

2. Include a similar parametric variation of mt when estimating uncertainties on the 
interpretation of any search in the region close to the EWSB boundary, since it can 
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change the kinematics of SUSY cascade chain decays, and therefore the efficiencies. It 
is intended that the theoretical uncertainty will be reduced in the future by including 
higher order effects from Refs. [32, 33] in SQFTSUSY. 

3. Use a recent version of SOFTSUSY. In version 3.3.6, we have improved the numerical 
convergence of SOFTSUSY near the boundary of EWSB. We have also improved the 
diagnosis of an accurate solution of the RGEs and boundary conditions by using the 
difference between the value of the Z mass predicted and the experimental value. 

We hope that ATLAS and CMS, as well as others, will adopt these recommendations as 
well as the Les Houches Recommendations [39], leading to more robust interpretations of 
the searches. 
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